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I. INTRODUCTION 

Very recently, the Belle Collaboration announced two 
charged bottomonium-like states Zi,(10610) and Z/,(10650). 
These two states were observed in the invariant mass spectra 
of htinP)^ (n = 1,2) and TimS)^ (m = 1,2,3) of the cor- 
responding T(55) -> hdnP)n + n- and T(55) -> f(niS)n + n- 
hidden-bottom decays |ll]. With the above five hidden-bottom 
decay channels, Belle extracted the Zi,(10610) and Z/,(10650) 
parameters. The obtained averages over all five channels are 
M z „(io6io) = 10608.4 + 2.0 MeV/c 2 , r Zi (io6io) = 15.6 + 2.5 
MeV/c 2 , M Zb(m50) = 10653.2 + 1.5 MeV/c 2 , r Zb(m50} = 
14.4 + 3.2 MeV/c 2 Q|]. In addition, the analysis of the angu- 
lar distribution indicates both Z;,(10610) and Z/,(10650) favor 
F(JP) = 1+(1+). 

If Zfe(10610) and Zb( 10650) arise from the resonance 
structures, they are good candidates of non-conventional 
bottomonium-like states. The masses of the J PC — 1 ++ and 
J PC = l + ~ bbqq tetraquark states were found to be around 
10.1 ~ 10.2 GeV in the framework of QCD sum rule for- 
malism |2], which are significantly lower than these two 
charged Z\, states. Therefore, it's hard to accommodate them 
as tetraquarks. If comparing the experimental measurement 
with the BB* and B*B* thresholds, one notices that Z/,(10610) 
and Z/,(10650) are close to thresholds of BB* and B*B*, re- 
spectively. One plausible explanation is that both Z/,(10610) 
and Zi,( 10650) are either B*B* or B*B* molecular states re- 
spectively. 

Before the observations of two charged Z/,(10610) and 
Zi,( 10650) states, there have been many theoretical works 
which focused on the molecular systems composed of B (t> 
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and meson pair and indicated that there probably exist 
loosely bound S-wave B*B* or B*B* molecular states |3j,|4[]. 
To some extent, such studies were stimulated by a series of 
near-threshold charomonium-like X, Y, Z states in the past 
eight years. 

Molecular states involving charmed quarks were first pro- 
posed by Voloshin and Okun more than thirty years ago Jit]. 
Later, De Rujula, Georgi and Glashow speculated 1^(4040) as 
a D*D* molecular charmonium [6]. Tonqvist calculated the 
possible deuteron-like two-meson bound states such as DD* 
and D*D* using the quark-pion interaction model J7|,|8t]. The 
observations of X(3872), three charged charomonium-like 
states Z + (4350), Z+(4050), Z + (4250) and r(4140), 7(4274) 
etc. again inspired theorists' interest in the molecular system 
composed of charmed meson pair (see Refs. l3l H l9l- [37ll for 
details). 

As the first observed charged bottomonium-like states, 
Z;,(10610) and Z fc (10650) have attracted the attention of 
many theoretical groups. The authors discussed the spe- 
cial decay behavior of the / = 1 S-wave BB* and B*B* 
molecular states based on the heavy quark symmetry in 
Ref. JH. Chen, Liu and Zhu found that the inter- 
mediate Z fe (10610) and Z fc (10650) contribution to T(5S) -> 
T(2S )n + 7i~ naturally explains Belle's previous observation of 
the anomalous Y(25 )n + n~ production near the peak of Y(55 ) 
at yfs = 10.87 GeV fH, where the resulting dr(T(5S) -> 
T(2S)n + 7T-)/dm„*„- and dF(T(5S) r(2S)n + 7T-)/dcos0 
distributions agree with Belle's measurement after inclusion 
of these Z/, states l39ll . The authors of Ref. ll4lll tried to repro- 
duce the masses of Z/,(10610) and Z/,(10650) using a molec- 
ular bottomonium-like current in the QCD sum rule calcu- 
lation. Yang et al. studied the mass spectra of the S-wave 
[bq][bq], [bq]*[bq], [bq]*[bq]* in the chiral quark model and 
indicated that Z/,(10610) and Z/,(10650) are good candidates 
of the S-wave BB* and B*B* bound states Bugg pro- 

posed a non-exotic explanation of Zt(10610) and Z;,(10650), 
which are interpreted as the orthogonal linear combinations of 
the qq and meson-meson states, namely bb+BB* and bb+B*B* 
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114 311 . respectively. Nieves and Valderrama suggested the pos- 
sible existence of two positive C-parity isoscalar states: a 
3 Si - 3 D\ state with a binding energy of 90-100 MeV and 
a 3 P state located about 20-30 MeV below the BB* threshold 
iPPill . Unfortunately, the quantum number of the above states 
does not match those of these two charged Zj states. Danilkin, 
Orlovsky and Simonov studied the interaction between a light 
hadron and heavy quarkonium through the transition to a pair 
of intermediate heavy mesons. Based on the above coupled- 
channel effect, the authors discussed the resonance structures 



for the isovector states, and 



Using the chromomagnetic 
461 1 discussed the possibility 



close to the B ( * ] B* threshold Q45| 
interaction, the authors of Ref. 
of Zi(10610) and Zfe(10650) being tetraquark states. In con- 
trast, the bbqq tetraquark states were predicted to be around 
10.2 ~ 10.3 GeV using the color-magnetic interaction with 
the flavor symmetry breaking corrections j47il . consistent with 
the values extracted from the QCD sum rule approach 01 . 

As emphasized in Ref. J39ll . future dynamical study of the 
mass and decay pattern of the S-wave BB* and B*B* molec- 
ular states are very desirable. In this work, we perform more 
thorough study of the BB* and B*B* systems using the One- 
Boson-Exchange (OB E) model. Different from our former 
work in Refs. ylQl, we not only consider S-wave interaction 
but also include D-wave contribution between B w and B ( *K 
Such a study will be helpful to answer whether the BB* and 
B*B* molecular bottomonia exist or not. 

This paper is organized as follows. After the introduction, 
we present the formalism of the study of the BB* and B*B* 
systems, which includes the relevant effective Lagrangian and 
coupling constants, the derivation of the effective potential of 
the BB* and B*B* system, the numerical results etc.. Finally, 
the paper ends with the discussion and conclusion. 



II. DEDUCTION OF EFFECTIVE POTENTIAL 

A. Flavor wave function 

We list the flavor wave functions of the BB* and B*B* sys- 
tems constructed in Refs. |H13]. The BB* systems can be 
categorized as the isovector and isoscalar states with the cor- 
responding flavor wave functions 



|Z^,[J]°) 



= 

V2 V 



"> + \B*°B* )) (4) 



for the isoscalar state. In the above expressions, the super- 
scripts T and S in Eqs. (H}-© are applied to distinguish the 
isovector and isoscalar states, respectively. The total angular 
momentum of the S-wave B*B* systems is / = 0, 1,2. Thus, 
we use the extra notation [J] in Eqs. ©-dill to distinguish the 
B*B* systems with different total angular momentum J. 

Belle indicated that both Z fc (10610) and Z b ( 10650) belong 
to the isotriplet states. If Z fc (10610) and Z/,(10650) are the BB* 
or B*B* molecular states respectively, they should correspond 
to Z^ t and Z^ t [l] in Eqs. ([TJ and (f3j, respectively. Since 
Z/,(10610)° is of C-odd parity, i.e., C = -1, thus the coefficient 
c — +1 is taken in Eq. (fTJ. The choice of the coefficient c — -1 
and C = +1 leads to X(3S72) and its partners, where X(3872) 
corresponds to Z^, listed in Table. U 

In Table|I] we summarize the quantum numbers of the states 
when we discuss whether there exist the BB* and B*B* molec- 
ular states. Moreover, we extend the same formalism to study 
the DD* and D*D* systems, where the flavor wave function of 
the DD* and D*D* systems can be obtained with replacement 
_> S« an d jj« _> £>(*). 



TABLE I: A summary of the BB*, B*B*, DD*, D*D* systems. If tak- 
ing c = - 1 in Eqs. (0 and l[2j, we obtain the flavor wave functions of 
Z nE i and Z„ s , , which are the partners of Z nS ' and Z nS , respectively. 



BB*/B*B* systems DD*/D*D* systems 


l G (J K ) 


^BB' 


7 (T) 
DD* 


l + d + ) 


^BB' 


7 (S) 
^DD* 


<ro + -) 






L(0 + ), L(2 + ), 1 + (1 + ) 




Z%,U] 


+ (0 ++ ),0 + (2 ++ ),0-(l + -) 


7 (T) ' 
BB* 


7 {T) ' 
^DD* 


L(l + ) 


7 (S) ' 
BB* 


7 (S) ' 
DD* 


+ (l ++ ) 



P%T> = ^B'+B^+cB+B* ), 
\^l:)-U\B*-~B") + cB-B*% (1 ) 
(\B* + B-) - B*°B°) + c(B + B*- - B°B*°) , 



(\B* + B-) + B*°B°) + c(B + B*- + B°B* Q ) 



(2) 



where c = ± corresponds to C-parity C = + respectively [3, 
4]. The flavor wave functions of the B*B* systems can be 
constructed as 



|z£U j n = \B* + B*") 
|Z^,[J]-) = \B*-B*°) 



(3) 



\Z { ll,V] ) = ^_{\B* + B*-)-\B* B* Q )) 



B. Effective Lagrangian and coupling constant 

In order to obtain the effective potential of the BB* and 
B*B* system, we employ the OBE model, which is an ef- 
fective framework to describe the BB* or B*B* interaction by 
exchanging the light pseudoscalar, scalar and vector mesons. 
In terms of heavy quark limit and chiral symmetry, the inter- 
actions of light pesudoscalar, vector and scalar mesons inter- 
acting with S-wave heavy flavor mesons were constructed in 
Refs. JMSF^ 



-Chhp = ig{H™y,A* ba y 5 l 

+ ig{Hf^ ah y 5 Hf\ 



(5) 
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+iA(Hf ) a, v F> iV {p)Hf>) 
tun, = gs (H^crH^) +gs (HfcrH^), 



(6) 
(7) 



where the multiplet field H (& is composed of the pseudoscalar 
P and vector P* with P ( * )T = (D w+ , D w0 ) or (B ( * )0 , fi w ~). 
And H (e) and are defined by 



[p^y-nrs], 



(8) 
(9) 



Here, H = y Wy and v = (1, 0). 

As given in Refs. I28ll54ll . the anti-charmed or bottom me- 
son fields P { * )T = (D (t) -, 5 ( * )0 ) or (fi ( * )0 , fi w+ ) satisfy 



(10) 



The multiplet field with the heavy antiquark can be de- 
fined as 



= C(CHrc-y C 1 = [P7r, - Pays] 



l-fi 



,(11) 



(G) 



[^+^75]. 



(12) 



If considering the following charge conjugation transforma- 
tion, 



c-fcr 1 = f, c%<r' = -<vl 



= CffC- 1 = -ff T , 



(13) 



one obtains the Lagrangian relevant to the mesons with heavy 
antiquark Q which is converted from the one related to the 
meson with heavy quark Q, where the Lagrangians are given 
in Eqs. ®-(7} jfHH. In the above expressions, the P(P) 
and p*(P*) satisfy the normalization relations {Q\P \Qq(0~)) = 

<0|^|<2<7(0-)> = v^and<oip;ie^(i-)> = «#J|fig(l-)> = 

e A , y/Afc. The axial current is A* = \{^d^ - %d^) = 



+ ■■■ with £, = exp(i¥/f„) and f„ = 132 MeV. 

PL = igv%J V2, />(p) = dppy - dyPf, + [p^, p v ], F'^ip) = 
- dypf, - Iff!, p v ] and g v = m p /f n , with g v = 5.8. Here, 
P and V are two by two pseudoscalar and vector matrices 



vr 



_ i n + jl 

V5 V5 / 



V2 V2 



p 
V3 



(14) 



(15) 



V2 ^ 



By expanding Eqs. (O-Q), one further obtains the effective 
Lagrangian of the light pseudoscalar mesons P with the heavy 



flavor mesons 

■L<p"p-i 



Jn 

+i^s a , vA v a PT t Pl A d v P ab , 

J" 



L<F» 



r? 



Jn 



(16) 



(17) 



The effective Lagrangian depicting the coupling of the light 
vector mesons V and heavy flavor mesons reads as 



(18) 



-Cppv = - ^l2pg v P b Plv ■ Y ba + y/2/3g v PlP b v ■ Yob, 

= -2y/2Ag v v A s A , afi (P b PT f +P2'P^(d a ¥' 3 ) ba 
-2 yl2Ag v v A e A , al3 {pfp h + PffiXd^U, 

(19) 

-Cp-p-y = ^2(3gvPl-P^v-Y ba 

-il ^llAgvP^pfid^v - dyVJia 

- y/2/3g v P*JP* b v ■ Yob 

-a yf2A gv p*fp* b v {d^ v - dyY^u. (20) 

The effective Lagrangian of the scalar meson cr interacting 
with the heavy flavor mesons can be expressed as 

Lv9a = -2gP b P\<T-2gP b Pl<T, (21) 

£p, P , a = 2g s Pl-P*J<r + 2g s Pl-P*Jo: (22) 

As shown in Eqs. (fT6b-(f20l>. the terms for the interactions 
between the anti-heavy flavor mesons and light mesons can 
be obtained by taking the following replacements in the cor- 
responding terms for the interactions between the heavy flavor 
mesons and light mesons: 



1 — > -v, a 

p: -> Pt 



P*' 



p;, p 



b, b — > a, 
P -> ~P\ 

-P. 



g = 0.59 is extracted from the experimental width of D* + 
115511 . The parameter /S relevant to the vector meson can be 
fixed as (3 — 0.9 by the vector meson dominance mechanism 
while A = 0.56 GeV~' was obtained by comparing the form 
factor calculated by light cone sum rule with the one obtained 
by lattice QCD. As the coupling constant related to the scalar 
meson cr, g, = g n /(2 V6) with g„ = 3.73 was given in Refs. 



C. Effective potential 

With the above preparation, we deduce the effective poten- 
tials of the BB* and B*B* systems in the following . Generally, 
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the scattering amplitude iAi(J, Jz) is related to the interaction 
potential in the momentum space in terms of the Breit approx- 
imation 



vn,-2M ; n/2M 



where M, and My denote the masses of the initial and final 
states respectively. The potential in the coordinate space 'V(r) 
is obtained after performing the Fourier transformation 



-vf B< *'(r) 



= J ^' p ^'*W^V,4x (23) 

where we need to introduce the monopole form factor (FF) 
J r (q 2 ,m^ E ) = (A 2 - m|)/(A 2 - q 2 ) to reflect the structure ef- 
fect of the vertex of the heavy mesons interacting with the 
light mesons, mg denotes the exchange meson mass. For 
q 2 — > we can treat FF as a constant while for A » m FF 
approaches unity. The behavior of FF indicates JH (1) when 
the distance becomes infinitely large, the interaction vertex 
looks like a perfect point corresponding of the constant FF; 
(2) when the distance is very small, the inner structure would 
manifest itself. In reality, the phenomenological cutoff A is 
around one to several GeV, which also plays the role of regu- 
lating the effective potential. 

In this work, we consider both S-wave and D-wave interac- 
tions between and mesons. In general, the BB* and 
B*B* states can be expressed as 



|flfl*( 3 Si) 
|fiB*( 3 £>i) 



. I44.ro) = 



B*B*( 3 Si)j ' 

B*B*( 3 D 1 )j 

B*B*( 5 D\)\ 



,|Z^.[2]) = 



B*B*( l So)) 


\ 


B*B*( 5 D )j 


/ 


B*B*( 5 S 2 ) 


>] 


B*B*{ l D 2 ) 


> 


B*B*( 3 D 2 ) 


> 


B*B*( 5 D 2 ) 


>) 


(24) 



with a = S,T, where we use the notation 2s+l Lj to denote the 
total spin S , angular momentum L, total angular momentum 
/ of the BB* or B*B* system. Indices S and D indicate that 
the couplings between B* and B* occur via the S -wave and 
D-wave interactions, respectively. 

Thus, the total effective potentials of the BB* and B*B* sys- 
tems are 



v 4"i.[ J i 

V Total 



z (a) (') 
BB* 



")| £ <Vf> (r)|zg. W ), (25) 

E=n,r},cr,p,(jj 

4%.[ J ])| Z V™(r)\Z%m), (26) 

E=Tr,T],o~,p,bj 



which are 2x2 and (J + 2) x (J + 2) matrices respectively. We 
impose the following constraint 



\bb*( 



25 + 1 



Z r JM 
bn,L 



(27) 



\b*b*( Lj 



c 



Sm s ,Lm L ^ \m,lm ,(z ri e " 2z -'"£ 



Sm s jri m v 



m,m ,mL,tns 



(28) 



to the effective potential obtained from the scattering ampli- 
tude. Ci M T , Cl M , and Cf'"i , are the Clebsch-Gordan 

\m,Lmi? Sms,Lmi lm,lm' 

coefficients. Yi mL is the spherical harmonics function. The 
polarization vector for the vector heavy flavor meson is de- 
fined as e 1 " = +-^=(e"' ± ie'") and e™ = ef . Here, the polariza- 
tion vector in Eqs. d27l>-(f28b is just the one appearing in the 
effective potentials which will be presented later. 



1. The BB* system 

The general expressions of the total effective potentials of 
the isoscalar and isovector BB* systems are 



1 



j_ y Direct _|_ j_ yuirect _j_ j_ j _ 



l 



/Cross 



_j__yCross _ 2^/Corss _j_ ^yCross 



3 1 r/ 

T/Direct . _ T/Direct . _wDirect . __| /-T/Cross 

* <r 2 2 ~ 



+ _yCross + 6y Coi SS + 2y Cor SS 



(29) 



(30) 



where the subpotentials from the n, rj, cr, p and u meson ex- 
changes are written as 



V 



Cross 



r- 

.In 



l 



-(e 2 • el)Z(A 2 ,m 2 ,r) 



+ -S (r, e 2 , e\)T(A 2 , m 2 , r) 



Cross 



f 2 

■In 



-(e 2 ■ el)Z(A 3 ,m 3 ,r) 



V 1 



V, 



Direct 



+ -S(?,e 2 ,el)T(A 3 ,m 3 ,r)], 
-8 2 s(^2 ■ <)T(A, m a , r\ 
~fg 2 v(e 2 -el)Y(A,m p ,rX 
2 , 

-(e 2 ■ e 3 )Z(A , m , r) 



2A 2 gl 



- -S (r, e 2 , 6pr(A , mo, r) 



V, 



Direct 



--B 2 g\(e 2 - e\)Y{A,m^r\ 



/Cross 



= 2A 2 g 2 y 



1 



'-(e 2 ■ e')Z(Ai,rai,r) 



--5(r, e 2 , 63)r(Ai,mi,r) 



(31) 



(32) 
(33) 
(34) 



(35) 
(36) 

(37) 



In the above expressions, we define 

A\ = A 2 - {m B > -m B ) 2 , m\=m\- (m B * - m B ) 2 , 
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A 

A 2 



3 = A - im B > - m B ) 

2 


2 
1 



= m 2 - (m B , - m B f, 



(m B > 
{m B > 



■ m B f, 

■ m B f, 



(m B . - m B ) , 
- (m B > - m B ) 2 , 



and S (r, a, b) = 3(r ■ a)(r • b) - a b. Additionally, functions 
F(A, m, r), Z(A, m, r) and T(A, m, r) are defined as 



Y(A,m E ,r) = 
Z(A,m E ,r) = 
T(A,m E ,r) = 



1 

4nr 



. A 1 - mi 
- <T Ar ) - 



8ttA 



(38) 



1 d 



v'Y(A,m E ,r) = -—^ — Y{A,m E , r \39) 
r L or or 



d 1 d 
r^--^-Y(A,m E ,r). 
or r or 



(40) 

(S ) 



In Eqs. ([29])-([30]), c = +1 corresponds to the Z*V and 
states including these two charged Zj states observed by Belle 



collaboration while taking c = -1 corresponds to the Z 

7 f0 ' 



<T) ' 
BB 



and Zl^ states which are partner states of X(3872). 

BB . . 

As indicated in Eq. ( 124l i. we consider both S-wave and D- 
wave interactions between the B and B* mesons. Finally the 
total effective potential can be obtained by making the replace- 
ment in the subpotentials 



(f2 ■ e\) 
(e 2 ■ e\) 



1- 


1 

1 













S(r, ei, el) 



- V2 
• V2 1 



which results in the total effective potential of the BB* system, 
i.e, a two by two matrix. 

The effective potential of the DD* system is similar to that 
of BB* system. The 77, <x, p and u> meson exchange potentials 
of DD* system can be easily obtained by replacing the param- 
eters for the BB* system with the ones for DD* system. Since 
the mass gap of m* D and mo is larger than the mass of n, which 
is different from the case of the BB* system, the n exchange 
potential of the DD* system is 0,0] 



V 



Cross 



r- 

.In 



l -ie 2 -e\)Z° D \A A ,m A ,r) 



+ -Sir,e 2 ,e< 3 )Tz" u (A 4 ,m 4 ,r) 



(41) 



where 



Y? D 'iA 4 ,m 4 ,r) 



1 

Anr 



,-A4r 



r(A 2 + m 2 ) _ 
2A 4 



+ cos(razir) 



(42) 



Z° D \A A ,m A ,r) = v 2 yf D *(A 4 ,m 4 ,r) = 4|"'- 2 |- 

r L or or 

xy™'(A 4 ,m 4 ,r), (43) 

Q \ Q 

T° D 'iA A ,m A ,r) = r Y ° D ' iA A ,n M ,r). (44) 

or r or 

In the present case, the parameters A 4 and m 4 are defined as 
A 4 = VA 2 " im D - - m D ) 2 , (45) 
m 4 = ^im D , -m D ) 2 -ml. (46) 



2. The B*B* system 

For the isoscalar and isovector B*B* systems, the general 
expressions of the total effective potentials are 

<V Z S* M = W a - l -W p + l -W^- l -W n + l -W ll , (47) 
111b 

W <T + \w p + \v a + \w x +\w n , (48) 
2 2 1 6 



<V Z B-Z> 



respectively, where the n, i], cr, p and a> meson exchanges 
can contribute to the effective potentials. The corresponding 
subpotentials are expressed as 



W„ = - 



r- 

J IT 



1 , , 

-(ei x e') • (e 2 x el)Z(A, m K , r) 



--5(r,ei x e[,e 2 x ehT(A,m„,r) 



1 



-(£1 x e\) ■ (e 2 x c;)Z(A, nir,, r) 



fn 

-is(r,ei xel,e 2 xel)TiA,m, ! ,r) 
-g 2 s iei-el)ie2-ehYiA,ma.,r\ 



1 



(49) 

(50) 
(51) 



2J3 2 gliei-ehie 2 -ehYiA,m p ,r) 
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1 



-ie\ x e\) ■ ie 2 x eJ)Z(A, m p , r) 



W„, 



--Sir,ei x e' v e 2 x e\)TiA,m p , r) 
= —hfgliex ■ e\)ie 2 • eJ)7(A, m w , r) 



(52) 



-M 2 gl 



-iei x • (e 2 x e^)Z(A, m w , r) 



--Si?,ei x e^,e 2 x eJ)r(A, m^, r) 



(53) 



Here, the definitions of F(A, m, r), Z(A, m, r), T(A, m, r) and 
5 (r, a, b) are given in Sec. Ill C 1 1 

In this work, we consider both S-wave and D-wave interac- 
tions between the B* and B* mesons, which are illustrated in 
Eq. d24b . Thus, the total effective potential of the B*B* with 
J — 0, 1, 2 is 2x2, 3x3, 4x4 matrices, which can be obtained 
by replacing the corresponding terms in the subpotentials, i.e., 



(61 • el)ie 2 ■ el) 



(61 x e\) ■ ie 2 x e\) 



Sir,e\ x 63,62 x eh 



1 
1 



2 
-1 

V2 

V2 2 



(54) 



(55) 



(56) 
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for the B*B* states with 7 = 0, 



(e, ■ e\)(e 2 • e\) 



1 
1 
1 



(ei x el) ■ (e 2 x el) 



5(r,ei XeJ,e 2 XeJ) 



1 
1 
0-1, 

( -V2 0^ 
-V2 1 
1 , 



(57) 



(58) 



(59) 



for the ZTZT states with J = 1, and 



(e, ■ eh(e 2 • e+) 



(ej x el) ■ (e 2 x el) 



10 
10 
10 
1 J 

-1 ' 

2 

10 

0-1, 



(60) 



(61) 



S(r,£i x e^,e 2 x el) 



o 

[u 2 
I V 5 V7 



0-^= 

V7 



0-10 







(62) 



for the B*B* states with J — 2. 

The potentials of the D*D* system and B*B* system have 
the same form. We only need to replace the parameters for 
the B*B* system with the ones for the D*D* system. 



III. NUMERICAL RESULTS 

With the obtained effective potentials, we can find 
the bound state solution by solving the coupled-channel 
Schrodinger equation. Corresponding to the systems in Eqs. 
(E5)-(ra, the kinetic terms for the Z ( "\ ' and Z ( ' v) - [J] (J = 

DD* D* B* 



0, 1 , 2) systems are 



K {a] > = diag 



A 
2tni ' 



A 2 
2fh i 



(63) 



diag 
diag| 
diag 



A 


A 2 ' 


2m 2 ' 


2m 2 1 


A 


A 2 


2m 2 ' 


2m 2 ' 


A 


A 2 


2ihi ' 


2fhi ' 



A 2 



At 



At 



(64) 
(65) 
{66) 



respectively. Here, A = ji§j.r 2 jp, A 2 = A - ^. mi = 
m gin b* / {m B + mB-) and fii2 = niB'/2 are the reduced masses of 
the zS and systems, where nis and ;7Zb» denote the masses 



of the pseudoscalar and vector bottom mesons 15911 . respec- 
tively. Of course, the kinematic terms for the DD* and D*D* 
systems are of the same forms as those for the BB* and B*B* 
systems, where we replace the mass of Z) w with that of B ( *\ 

In this work, the FESSDE program |[56l l57ll is adopted 
to produce the numerical values for the binding energy and 
the relevant root-mean-square r with the variation of the cut- 
off in the region of 0.8 < A < 5 GeV. Moreover, we also 
use MATSCE JH, a MATLAB package for solving coupled- 
channel Schrodinger equation, to perform an independent 
cross-check. 

Throughout this work, we will first present the numerical 
results of the obtained bound state solutions when all types of 
the one-meson-exchange (OME) potentials are included. The 
one-pion-exchange force contributes to the long range inter- 
action between the heavy meson pair, which is clear and well- 
known. In contrast, the scalar and vector meson exchanges 
are used to mimic the intermediate and short range interaction 
between the heavy mesons, which are not determined very 
precisely. In order to find out whether the existence of the 
possible bound molecular states is sensitive to the details of 
the short-range interaction, we will also study the case when 
only the one-pion-exchange (OPE) contribution is considered. 
In the following illustration, we use OME and OPE to distin- 
guish such two cases. If the OPE force alone is strong enough 
to form a loosely bound state, such a case is particularly inter- 
esting phenomenologically. 



A. The BB* and DD* systems 



In the following, we first present the numerical results for 
theZis andZi,2 states where zl 7 ^ corresponds to Zh( 106 10) 
observed by Belle HI]. As shown in Table |H there exist two 
systems with c = — 1 and C — +1 in the flavor wave functions, 
where are marked as Zi2, and Z^, . 

DD DD 

1. In Table HD we present the numerical results of the ob- 
tained bound state solutions in both OME and OPE 
cases. We find the bound state solutions for the two 

• (S} (S} ' 

isoscalar Z„= and Z\l with reasonable A values (A ~ 

DD* DD* 



1 GeV), which indicates the existence of the z' 5 -' 

DD* 



and 



■'BB* 



molecular states. 



2. For the Zr„l, state, we also find the bound state so- 

BB 

lution with A around 2.2 GeV. Our result shows that 
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Z could be as a molecular state with a very shallow 

DD 

binding energy. In addition, its binding energy is not 
strongly dependent on A. Thus, it is quite natural to in- 
terpret Z/,(10610) as a BB* molecular state with isospin 
7=1. 

(7") ' 

3. For the Z - system, the bound state solution can be 
found in the region A > 4.7 GeV. To some extent, the 
value of A for the Z^ 3 ? seems a little large compared to 
1 GeV. 

4. We also discuss the case when we only consider the 
OPE potential. For the {Z^,, Z ( n r > ', Z ( n s ' '} or Z ( j2„ we 

DD DD DD DD 

need to decrease or increase the A value to obtain the 
same binding energy as that from OME. The one pion 
meson exchange potential indeed plays the crucial role 
in the formation of the BB* bound states. 



TABLE II: The obtained bound state solutions (binding energy E and 
root-mean-square radius trms) for the BB* systems. Here, we discuss 
two situations, i.e., including all one meson exchange (OME) contri- 
bution and only considering one pion exchange (OPE) potential. 





OME 


OPE 


I a (J pc ) 


State A E (MeV) trms (fm) 


A E (MeV) trms (fm) 


1 + (1 + ) 


2.1 -0.22 3.05 
2.3 -1.64 1.31 
2.5 -4.74 0.84 


2.2 -8.69 0.62 
2.4 -20.29 0.47 
2.6 -38.54 0.36 


r(i + ) 


4.9 -0.14 3.64 
Z { pj 5.0 -0.41 2.45 
5.1 -0.85 1.80 


4.5 -17.79 0.56 

4.6 -22.65 0.52 

4.7 -28.29 0.48 


o-(i + -) 


1.0 -0.28 3.35 
Zfl 1.05 -1.81 1.71 

1.1 -5.36 1.18 


1.8 -10.09 0.96 

1.9 -15.11 0.84 
2.0 -21.53 0.76 


+ (l ++ ) 


0.8 -0.95 1.84 
Zfl ' 0.9 -6.81 0.91 
1.0 -19.92 0.65 


1.0 -7.68 0.82 

1.1 -15.30 0.65 

1.2 -26.53 0.53 



We extend the formalism in Sec. HI] to study the DD* sys- 
tems. As shown in Table HTT1 we can exclude the existence of 
the Z^, and Z^, since we do not find any bound state solu- 
tion for the Z^ and Z^, states. For the two isoscalar zj^, 

(S ) ' 

and ZV:, , there exist loosely bound states with reasonable A 
values. If only considering the OPE exchange potential, we 
notice: (1) the bound state solution of the Z^ appears when 
A ~ 4.6 GeV, which largely deviates from 1 GeV; (2) there 
still does not exist any bound state solution for 7^ D L ; (3) for 

dd' anc ^ Db' ' becomes larger in order to find the bound 
state solution. The comparison between the OME and OPE 
results also reflects the importance of one pion exchange in 



the DD* systems. We need to specify that Z^l' with + (l ++ ) 
directly corresponds to the observed X(3872) 16011 . 

The BaBar Collaboration measured the radiative de- 
cay of X(3872) and found a ratio of B(X(3872) 
i//(2S )y)/fi(X(3872) -> J/if/y) = 3.4+1.4 JUl], which contra- 
dicts the prediction with a purely DD* molecular assignment 
to X(3872) (12]. However, very recently Belle reported a new 
measurement of the radiative decay of X(3872), where only 
the decay mode X(3872) — > J /(fry was observed and the upper 
limit fi(X(3872) f (2S)y)/B(X(3S72) J I iffy) < 2.1 was 
given 16211 . The inconsistence between the Belle and BaBar 
results indicate that the study of X(3872) is still an important 
research topic. Our numerical results suggest that the mass 
of the loosely bound molecular state Z^ t is consistent with 
that of X(3872). The assignment of X(3872) as a molecular 
candidate is still very attractive. 

TABLE III: The obtained bound state solutions (binding energy E 
and root-mean-square radius ?rms) for DD* systems. Here, we dis- 
cuss two situations, i.e., including all one meson exchange (OME) 
contribution and only considering one pion exchange (OPE) poten- 
tial to BB* systems. 





OME 


OPE 


I a (J FC ) 


State A E (MeV) trms (fm) 


A E (MeV) trms (fm) 


1 + (1 + -) 


7 (T) 
^DD' 


4.6 -0.85 1.46 

4.7 -3.42 1.17 

4.8 -7.18 0.93 

4.9 -12.40 0.75 


r(i ++ ) 


7 IT) 1 
^DD' - 




o-(i + -) 


1.3 

m 1.4 -1.56 1.61 

z 

DD * 1.5 -12.95 0.98 
1.6 -35.73 0.69 


3.4 -0.11 1.74 

3.5 -2.03 1.50 

3.6 -4.79 1.26 

3.7 -9.62 1.06 


+ (l ++ ) 


1.1 -0.61 
.„ , 1.2 -4.42 1.38 

Z 

DD " 1.3 -11.78 1.05 
1.4 -21.88 0.86 


1.7 -3.01 1.37 

1.8 -7.41 1.06 

1.9 -14.15 0.84 
2 -23.82 0.68 



B. The B*B* and D*D* systems 

The numerical results of the B*B* systems are presented in 
Table IIV1 which include the obtained binding energy and the 
corresponding root-mean-square radius. We find the bound 
state solution for all the B*B* states with reasonable A values: 

( T) 

1. A loosely bound state exists for Z B ,^[1] corresponding 
to the observed Z h ( 10650) with A slightly above 2 GeV. 
Only considering the OPE potential, the obtained bind- 
ing energy becomes deeper with the same A value. 
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2. In addition, the B*B* can form loosely bound molecu- 
lar states Zj" [0], Z%>[0], Z^Jl] and Z^[2] with 
very reasonable A values. Comparing the results be- 
tween OME and OPE cases, one notices again that the 
one pion exchange indeed is very important to form the 
B"B* bound state. 

3. For the Zj^ [2] state, the existence of the loosely bound 
state requires the value of A around 4.4 GeV. 



TABLE IV: The obtained bound state solutions (binding energy E 
and root-mean-square radius r RMS ) for B*B* systems. Here, we dis- 
cuss two situations, i.e., including all one meson exchange (OME) 
contribution and only considering one pion exchange (OPE) poten- 
tial to Bfi'systems. 





OME 


OPE 


I a (J pc ) 


State A E (MeV) rRMS (fm) 


A E (MeV) trms (fm) 


l + (0 + ) 


1.2 

1.4 -1.44 1.24 
1.6 -6.16 0.77 
1.8 -15.15 0.54 


1 

1.2 -0.32 1.53 
1.4 -5.69 0.78 
1.6 -18.82 0.50 


o-(o + -) 


0.9 

1 -0.81 2.11 

1.1 -9.98 1.02 

1.2 -35.16 0.70 


1 

1.2 -0.52 2.76 
1.4 -5.74 1.12 
1.6 -20.92 0.77 


1 + (1 + ) 


2.2 -0.81 1.38 
2.4 -3.31 0.95 
2.6 -7.80 0.68 
2.8 -14.94 0.52 


2 -2.17 1.15 
2.2 -8.01 0.68 
2.4 -19.00 0.48 
2.6 -36.36 0.38 


o-(r-) 


1. -0.01 2.07 

1.1 -5.50 1.17 

1.2 -21.76 -0.75 

1.3 -53.68 0.55 


1.4 -0.51 1.90 
1.6 3.65 -1.32 
1.8 -10.26 0.96 
2.0 -21.81 0.75 


l + (2 + ) 


4.4 -0.44 1.59 
4.6 -1.59 1.28 

Z<? fl .[2] 

4.8 -3.42 1.01 
5. -6.16 0.81 


3.6 -2.82 1.12 
3.8 -6.21 0.85 
4.0 -11.41 0.68 
4.2 -18.77 0.57 


0-(2 + -) 


0.8 -2.33 1.32 
2% ] B ,m 0.9 -10.45 0.84 
1.0 -27.14 0.63 


0.8 -1.81 1.48 
0.9 -5.64 1.01 
1.0 -12.28 0.76 



In the following, we also present the numerical results for 
the D*D* systems in Table M Our calculation indicates: 

1 . We find the bound state solutions for the three isoscalar 
states Z^JO], Z^Jl] and Z ( ^ D ,[2], where the cor- 



responding A is around 1 GeV. If only considering the 
OPE contribution for the Z^,[0], Z^Jl] states, we 
need to largely increase A value in order to obtain a 

(S) (S) 

loosely bound state. Here, either Zr D ,=, [0] or Zr = \2] 
could correspond to the observed 7(3930) by Belle J63I1 
and BaBar |640 , which is consistent with the conclusion 



inRef. OJ. 



2. There does not exist the bound state Z^^, [2] . The value 
of A is about 3.6 GeV in order to form a bound state 
Z^? 5 ,[0]. In the range 0.8 < A < 5 GeV, we cannot 

find the bound state solution for Z^,[l] in the OME 

case. Thus, we exclude the existence of the Z^^Jl] 
molecular state. 



TABLE V: The obtained bound state solutions (binding energy E and 
root-mean-square radius trms) for D*D* systems. Here, we discuss 
two situations, i.e., including all one meson exchange (OME) con- 
tribution and only considering one pion exchange (OPE) potential to 
D*D* systems. 





OME 


OPE 


I G (J PC ) 


State A E (MeV) r RMS (fm) 


A E (MeV) r RMS (fm) 


l + (0 + ) 


3.6 -0.94 1.74 
3.8 -6.16 1.00 

Zgk[0] 

4 -16.44 0.66 
4.2 -33.23 0.49 


2.8 -2.03 1.47 

2.9 -6.10 1.00 
3 -12.51 0.74 

3.1 -21.56 0.59 


<r(o + -) 


1.4 -1.72 1.62 

1.5 -17.98 0.88 

4 5 .Uo] 

1.6 -54.60 0.47 


3 -5.70 1.24 

3.1 -12.15 0.96 

3.2 -21.83 0.78 


1 + (1 + ) 


z£Ui] - 


4.7 -6.96 0.94 

4.8 -12.29 0.73 

4.9 -19.36 0.60 
5 -28.31 0.51 


o-(i + ~) 


1.3 

1.4 -3.44 1.44 

2&Ul] 

1.5 -16.57 0.90 

1.6 -41.25 0.66 


3.6 -9.91 1.01 

3.7 -15.25 0.87 

3.8 -22.07 0.76 

3.9 -30.53 0.68 


l+(2 + ) 






0-(2 + -) 


1.1 -0.61 1.72 

1.2 -7.50 1.19 

Z£U2] 

1.3 -19.22 0.89 

1.4 -35.93 0.73 


1.6 -3.89 1.28 

1.7 -9.64 0.98 

1.8 -18.38 0.77 

1.9 -30.71 0.64 
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IV. SUMMARY 

Stimulated by the newly observed bottomonium-like states 
Z/,(10610) and Z/,(10650), we have carried out a systemati- 
cal study of the BB* and B*B* system using the one boson 
exchange model in our work. We have considered both the S- 
wave and D-wave interaction between the B ( * } and B* mesons, 
which results in the mixing of the S -wave and D-wave contri- 
bution as discussed in Sec. [TT] Our numerical results indicate 
that the Z/,(10610) and Zi,(10650) signals can be interpreted 
as the BB* and B*B* molecular states with I G (J P ) = 1 + (1 + ) 
respectively. 

As a byproduct, we also predict the existences of six other 
BB* and B*B* bound states (see Table IVII ) within the same 
framework. We want to stress that the long-range interac- 
tion between the heavy meson pair arises from the one-pion- 
exchange force, which is clearly known. This OPE force alone 
is strong enough to form the above loosely bound molec- 
ular states, which makes the present results quite model- 
independent and robust. 

The observation of these Z/,(10610) and Z/,(10650) states 
shows that the hidden-bottom decay are very important de- 
cay channels, which is characteristic and helpful to the search 
of the molecular bottomonium. After taking into account of 
the phase space 159l 165146711 and the conservation of quan- 
turn number, the Z^„ Zjf-V, Z^.[0], Z^.[0], Z^.[l] and 

(S ) 

ZL 5 J2] molecular states can decay into 

(T( 1 S )!/, T(2S )rj, h b { 1 P)?j, m { 1 S )a>), 
{r(lS)a>,Xbo(lP)ri,Xbi(lP)v,Xb2(lP)v}, 

[t(\S)a>, Xb i{\P)Ti, TlbaS)T}}, 



{xbo(lP)co, T(1S)jj, T(25)/7, rj b (lS)cj, h b (lP)r]}, 

{r(lS)o),Xbi(lP)ri,Xb2(mv,Vb(lS)v}, 

respectively. The above modes can be used in the future 
experimental search of the partner states of Z/,(10610) and 
Z,(10650). 

We also extend our formalism to study the molecular char- 
monia. The observed possible molecular charmonia are listed 
in Table[VT] The possible hidden-charm decay channels of the 
molecular states z£]„ Z^, [0], Z^[l] and Z^. [2] are 

{r} c (lS)co,J/il,aS)Ti}, 

[j/lfr(J,T] c (lS)Tl}, 

lT) c (lS)oj,J/i/r(lS)T,}, 
[j/iKlS)a>, /7 C (15)7,}, 

respectively. Due to the limit of phase space, the hidden- 

(S) ' 

charm decays for the other one Z ^, molecular state are 
J/ifs(lS) or r/ c (lS) plus multi-pions. 
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BB* 
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7 {T) 
^DD* 
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o-(i + -) 


7 (S) 
^BB* 


/ 




7 (S) 
DD* 


/ 




1-(1 + ) 


7 m ' 

BB* 


X 




7 iT) ' 
^DD* 


X 




+ (l ++ ) 


7 (S) 1 
^BB* 


/ 




7 (S) ' 

^Db* 


/ 


X(3872) [60J 


L(0 + ) 




✓ 




Z D*D*M 


X 




+ (0 ++ ) 


Z b*bM 


/ 






/ 


7(3930) [65-67J 


1 + (1 + ) 




■/ 


Z b ( 10650) 




X 




o-(i + -) 




V 




z D*D*m 


/ 




l-(2 + ) 


Z b T * b *M 


X 




Z%.[2] 


X 




+ (2 ++ ) 




/ 




Z%,[21 
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7(3940) [65-67J 
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